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, The production of the Standard Model (SM) Higgs boson (H) plus one jet is compared with 

. that of the lightest scalar Higgs boson [h") plus one jet and that of the pseudoscalar Higgs boson 

' i^") plus one jet. The latter particles belong to the Minimal Supersymmetric Model (MSSM). We 

^S) I include both top and bottom quark loops to lowest order in QCD and investigate the limits of small 

quark mass and infinite quark mass. We give results for both the CERN Large Hadron Collider 
^ ' (LHC) and the Fermilab Tevatron. 

PACS numbers: 13.85.-t, 14.80.Bn, 14.80.Cp 
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. I. INTRODUCTION 

The Higgs boson is the cornerstone of electroweak symmetry breaking in the Standard Model (SM). Particle 
physicists around the world have made the search for the Higgs boson the top priority in high energy experiments. 
' However, there are several different candidate models in the Higgs sector. The Minimal Supersymmetric Standard 
' Model (MSSM), which is a special case of the Two Higgs Doublet Model (2HDM), is of particular theoretical interest. 
H The Standard Model Higgs boson has been experimentally excluded by LEP searches for e+e~ — ^ ZH if its mass is 

lighter than approximately 114 GeV/c^P|. In the MSSM, the particle spectrum includes five physical Higgs bosons; a 
JL light and a heavy neutral scalar (/i", i?°), two charged scalars (H^), and a CP-odd pseudoscalar (A°). The mass of the 
lightest scalar in the MSSM is excluded from being lighter than 91 GeV/ c^ , while the mass of the pseudoscalar M^o is 
, experimentally excluded from being lighter than approximately 92 GeV/c^. The ratio between the vacuum expectation 
■ values (VEVs) of the two neutral Higgs bosons of the MSSM is defined as tan/3 — V2/V1. For mtop = 174.3 GeV/c^, 
0.5 < tan P < 2.4 has been excluded by the LEP Higgs searches. A different value of the top quark mass will lead to 
different exclusion bounds on tan 
^ The total cross-section for scalar Higgs production including massive quark loops has been calculated at next-to- 

leading order (NLO) inperturbative QCD The corresponding calculation for Higgs production in the MSSM 

can be found in Ref . |a| . In the Heavy Quark Effective Theory (HQET) 0, S j the top quark mass is assumed 
to be much heavier than the Higgs boson mass and all relevant energy scales. Assu ming the HQET total inclusive 
cross-sections have been calculated at NLO for scalar^ and pseudoscalar productionfTlL fl^ and also at NNLO for 
scalar[ll[ll|ll[li| and for pseudoscalar [ll[ll|ll| production, see also|il|23|. The use of the HQET significantly 
simplifies the computation of higher order QCD effects and has been shown to accurately reproduce the exact NLO 
rate at the LHC for — > iJ 3, 5] for a Higgs mass less than 1 TeV/c^ if the LO massive results are multiplied by the 
NLO K- factor obtained in the HQET. 

In this paper we concentrate on the Higgs plus one jet {gg g^, qg (7$, and qq — > g^) production processes 
since they are important for the experimental detection of the Higgs. Here $ represents either the SM Higgs, H, or 
the MSSM scalars, h° and iJ", or the MSSM pseudoscalar A°. The production of the SM Higgs plus one jet process 
has been calculated exactly at LO in |^ with the inclusion of heavy quark loops. The production rate in the 
MSSM for the lightest scalar plus one jet was recently calculated in LO including SUSY loop effects, which can be 
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TABLE I: Higgs-fermion couplings in the MSSM and the dependence of the matrix element-squared on the couplings. Mt and 
Mb represent contributions from top- and bottom- quark loops, respectively. The a parameter is the angle that diagonalizes 
the CP-even Higgs squared-mass matrix. 



significant for light SUSY squarks and gluinosj2^. The NLO QCD corrections to the Higgs plus one jet process have 
only been computed in the HQET, since the full virtual corrections would require the evaluation of massive two- loop 
integrals for a 2 2 reaction. The differential cross-section for the production of a scalar Higgs boson plus one jet 
in the HQET at NLO has been calculated previously by [13, l20l 124 . l25l \2dL l27l | and the integrated rate was shown to 
increase substantially from the lowest order rate. The pseudoscalar case has been presented in [2^ and in p9l |. 

We present the calculation of the Higgs plus one jet process where we include both top and bottom quark loops 
with the full quark mass dependence. This is done for the SM Higgs and for the lightest scalar and pseudoscalar 
Higgs bosons of the MSSM. The contributions of loops with bottom quarks can be important for large values of tan /3 
in the MSSM. We also address the region of validity of the HQET predictions for these reactions. 

In Section ^ the limit of the partonic matrix elements in the HQET and in the small quark mass limit are 
explored. In Section lTTll the Higgs plus jet matrix elements are given and our computational techniques are described. 
Section llVI summarizes our notation for the hadronic differential cross-sections. Section Ivl contains numerical results 
for differential cross-sections at the Tevatron and LHC, as well as integrated results with cuts in transverse momemtum, 
Pt, and rapidity, y. Analytic results for the matrix elements are given in two Appendices. 



II. PARTONIC PROCESSES - HEAVY QUARK EFFECTIVE THEORY 



In the limit where the top quark mass is much heavier than all the energy scales in the problem, only the top 
quark coupling to gluons is numerically significant and this limit provides a good approximation to Standard Model 
Higgs production matrix elements. The HQET limit for scalar Higgs production has been extensively studied in the 
literature. This limit is especially useful for deriving higher order QCD corrections since the massive top quark loops 
that couple the Higgs boson to gluons reduce to effective vertices. The Feynman rules can be derived from an effective 
Lagrangian density 11 H S H [13, ilj , 
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where 5// = 1 in the Standard Model and v = 246 GcV. Oh generates vertices which couple the Higgs boson 
to two, threej_and four gluons. In the large rritop limit, the coefficient Ch can be evaluated as a power series in 
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where ai^\fj,'^) is evaluated at the scale /i^ in a 5 flavor scheme. 

For comparison, we consider a pseudoscalar Higgs boson with a coupling to fermions given by, 
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FIG. 1: Sample Higgs plus one jet diagrams. Figures la,b are the gg diagrams, Ic is the qg channel, and Id is the qq channel. 
All quarks contribute to the loops. The crossed and charge conjugate diagrams are not shown. There are a total of 12 gg 
diagrams and 2 for each of the qg and qq sub-processes. 



In the large mtop limit, the interactions of the pseudoscalar with gluons can be found from the effective Lagrangian 

[THHHlIll 

where GJ^" is the gluon field strength tensor. The process independent coefficient functions are 

Ga,=-^, CM^Oial). (5) 

We consider = 1 and the examine the differences between differential cross-sections for the production of a SM 
scalar Higgs boson and a pseudoscalar Higgs boson with the couplings of Eq. |3| when the bosons are produced in 
association with a jet. 

It is also of interest to compare the production rates for a Higgs boson plus a jet in the MSSM. The effective 
Lagrangians in this case are found by making the replacements in Eqs. Qand^l 

9H -^4 

9A^ct\ (6) 

where cf and cf° are given in Table I. (We neglect contributions from SUSY particles such as the bottom squarks 
and gluinos, and therefore assume that the SUSY particle masses are much larger than mtop and m$. These genuine 
SUSY contributions can be important for light squark and gluino masses [l^.) When the bottom quark becomes 
important, the HQET breaks down as a reliable calculational tool. This occurs in the MSSM when tan/3 becomes 
large and the bottom quark couplings are enhanced. 



III. PARTONIC PROCESSES - FULL THEORY 



There are three channels associated with Higgs plus one jet production: gluon fusion, quark-gluon scattering, and 
quark-antiquark annihilation. Representative Feynman diagrams are shown in Fig.^ At the LHC where \/S — 14 TeV 
the gluon fusion and quark-gluon channels are the most important, with the quark-antiquark channel adding a 
negligible amount to the process. However all three channels are important at the Tevatron where ^/S = 1.96 TeV. 
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FIG. 2: The squared matrix elements, \M\'^, evaluated at s = 4m| and w = i for the three different channels, {gg, qg, qq), 
normalized to the squared matrix elements in the HQET for scalar and pseudoscalar (with gA = 1) Higgs plus jet production. 
We include only the top quark loops. The solid line is the scalar, whereas the dashed line is the pseudoscalar. 
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The calculation of the matrix elements was carried out in both n = 4 — 2e dimensions and 4-dimensions. The 75 in 
the pseudoscalar calculation was treated using the Akyeampong-Delbourgo prescription|3^|33,|33| for the 75-matrix. 
In this scheme the 75 is exchanged for a Levi-Civita tensor contracted with four 7-matrices. After the trace, the 
tensor loop integrals were reduced to scalar integrals using the usual Passarino-Veltman|39j reduction techniques. 

A. Gluon fusion (gg (7$) 

The gluon fusion channel is the most important channel at the LHC. The momentum distribution in this process 
is assigned with all momentum incoming, 

5K"''^)+.9(pr''^^) -g(-pr''^^)+$(-P5), (7) 

where fj,i are Lorentz indices and Ui are color indices. The Mandelstam variables used in the partonic system are 

s = (Pl+P2)^ ^=(pl+P5)^ ^t=(p2+P5)^ Q^^ml. (8) 
The matrix elements, including the gluon polarization vectors, can be written 

M^^ - e^'{pi)e^HP2)e^'{P3)MZ.,^^- (9) 

The Ward-Takahashi identities let us check the gauge invariance of the sub-process. In the gluon fusion case, these 
can be written as 

giving us a strong check on the algebraic results. Analytic results for the matrix element squared fo r qq — > gA'^ are 
given in Appendix A, see also Appendix C in 5], while those for gg gH can be found in Refs. |2lL |22| . 

B. Quark-antiquark annihilation {qq — > g<l?) 

For this sub-process, the momentum, color, and Lorentz structure was assigned as follows 

q{pi) + q{p2) - + ^{-P^y (11) 

The matrix elements satisfy 

M^^^e>^-{p^)M%, P^'M%=Q. (12) 

Analytic results fo r qq gA? are given in Appendix B, see also Appendix C in while those for qq — > gH can be 
found in Refs. [21I l22||. The results for quark-gluon scattering can be found by crossing. 

C. HQET Matrix Elements 

The 4-dimensional color-spin averaged matrix elements for Higgs boson plus one jet production in the mtop 00 
limit are presented here for completeness. These matrix elements obey the same crossing relations as the full matrix 
elements, 

\M{S,i,u)\lg^^^ = -\Miu,i,s)\l^^g^. (13) 
The matrix elements in the large mtop HQET limit can be written [sl. I211 123 . l29| . 

Ei-^iL^.* = ^*4(iv?-T) Jfu ^''^ 

1 p + p 



El-^l^.-.*--^*^-^ (15) 
VlyU|2 -A (^c^-1) i^' + i' 
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where, 

-Ij-fe) (4™,)s5 (17) 

and gn — 1 for the SM and 5$ is given in Eq. |Blfor the MSSM. The bar imphes a sum and average over colors and 
spins. The exact matrix elements squared as compared with the HQET matrix elements are shown in Fig. [21 for both 
the SM scalar, which are in excellent agreement with the plots in |22| . and for a pseudoscalar with = 1. In this 
plot, the mass of the Higgs was set to m$ = 100 GeV/c^ and the mass of the top quark was varied. In these plots, 
two thresholds can be observed. Each threshold occurs when an imaginary part of the matrix elements turns on or 
off. If we examine Eq. IB14I for qq gA^ we clearly see that the imaginary part contains the difference of two step 
functions 

9{§ - iml^) ~ 9iM%, - 4m?„p), (18) 

so the first threshold occurs at 2mtop — M^o and the second at 2mtop = V§- Since we choose s — 4A/^o for the plot 
this implies that these thresholds occur at rritop/m^ = 0.5 and 1 respectively. The imaginary part is finite between 
these cusps. Similar phenomena occur in the other reactions. However when the squared matrix elements contain 
several terms the onset of the imaginary parts is not always visible. The reactions qg —^ g$ do not have s channels 
so they only have cusps at mtop/m^ = 0.5. Finally the gg — + g$ channels show both cusps. Note that the reason the 
cusps do not appear exactly at 0.5 and 1 is due to our choice of points in mtop/™*- 

These ratios show that when the heavy quark becomes heavier than mtop ~ the HQET is a reasonable 

approximation to the matrix elements with a top loop only. In the MSSM, however, the usefulness of the HQET is 
limited to small values of tan /? where the bottom quark contribution can be neglected. 



D. Small Quark Mass Limit 



When the quark mass in the loop is much smaller than the Higgs mass and the energy scale, the small quark mass 
limit m/ ^ is relevant. This is the case for the bottom quark contribution in the large tan/3 limit of the MSSM. 
The matrix elements in this limit behave as 

\M\' ^mjlog\m}/^^'), (19) 
where fi ^ ruf. Exact expressions in the small quark mass limit are given in Appendix B. 



IV. OBSERVABLES 

Generically, we can write a 2 ^ 2 differential observable as 

^':ll^ = T^Ei-^l'' (20) 

dt du Idtt ^ 

where the bar implies a sum and average over colors and spins. To relate the hadronic differential distributions to the 
partonic differential distributions we need to perform a convolution with the parton distribution functions. 
The hadronic process can be written as 

HiiPi) + H^{P2) ^ j(-P3) + *(-P5) (21) 
where the j represents the gluon or the quark jet in the sub-process of interest. In the hadronic system, we can write 

S^{Pi+P2)\ T^{Pi+p^)\ U^{P2+p^f. (22) 
This translates into the partonic system (with momentum fractions xi and X2) as 

pi = xiPi, p2 = X2P2, (23) 

s^xiX2S, i^xi{T-Q^) + Q^, u = X2{U - Q^) + (24) 
_ -U _ -xi{T-Q^)-Q^ 

Xl,n.i„- 2;2,ni„- ^ ^ _ q2 ' i^^i 
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where = m\ . The hadronic variables can be written in terms of the transverse momentum and rapidity 



T = Q2 ^ pI + Q2 coshy + \^Jpl + Q2 sinhy (26) 



= Q2 _ V^^p2 ^ Q2 cosh y - x/^y'p? + g2 sinh y. (27) 
The hadronic differential cross-section is, 

^''^ = T.t -f ^/fnx..?)/fnx2,.pi^^. (28) 
dT dU ^ Ai,„i„ Jx2,^in ^2 ■' dt du 

Upon further integration we obtain the single differential pt and rapidity distributions with the kinematic limits, 

_1 S~Q' _1 fl±SQ\ 

^""2 V5 ' ^'"""^2^1-5q;' ^ ' 

where Sq = JI-^^^^. (30) 



V. NUMERICAL RESULTS 



We present our calculations for the CERN LHC with VS* = 14 TeV and the Fermilab Tevatron with ^/S ~ 1.96 TeV. 
We use the CTEQ6.1L parton distribution functions 40] with Ag*-* = 165.2 MeV and a one loop running coupling 
constant with as{Mz) = 0.1298. For the differential distributions, the full kinematic rapidity and pt are used and the 
factorization and renormalization scales are set equal to, 

Mr = Ai/ = \/Q^+pI (31) 

We use pole masses with mtop = 174.3 GeV/c^ and mbot — 4.5 GeV/c^. For the integrated cross-section we require 
the Pt of the Higgs and the jet to satisfy pt,min > 30 GeV/c in the rapidity region \y\ < 2.5 and replace pt by pt^min in 
Eg. 1311 for the renormalization and factorization scales. 



A. Standard Model 



The transverse momentum distributions of the SM Higgs boson for all the separate channels are shown in Fig. O 
for the LHC. For a SM Higgs boson with Mh = 120 GeV/c^, the cross-section for Higgs plus one jet is approximately 
12.3 pb when both the top and bottom quarks are included in the calculation. Although the bottom quark contribution 
alone is only 0.05 pb, the top-bottom interference lowers the cross-section by approximately 8.25% from 13.4 pb when 
only the top quark is included, see ^Jl- This lowering of the cross-section may be visible at the LHC. As shown in 
Fig. 21 the full theory and the HQET agree very well at small to moderate pt for both the scalarj^ and the 
pseudoscalar differential distributions. 



B. Minimal Supersymmetric Standard Model 



The MSSM is a special case of the 2HDM. In the MSSM, the up- and down-type quarks become massive from 
different Higgs doublets and the ratio of the two VEVs is parameterized by tan (3 = V2 /vi . As shown in Table 
up- and down-type quarks couple differently to the Higgs bosons of the MSSM. The a parameter is the angle that 
is introduced to diagonalize the mass eigenstates of the CP-even Higgs squared-mass matrix to obtain the physical 
states. The program HDECAYji^ was used to determine the mass of the lightest scalar and the a mixing parameter 
once the values of Alj^o and tan/? were chosen. The SUSY Higgs mixing parameter was set to ji — 300 GeV/c^, the 
gluino mass to /Lt2 = 200 GeV/c^, all the SUSY breaking masses to 1 TeV/c^, and the soft breaking term to 1.5 TeV/c^. 

At the Tevatron, there is a very small signal for the SM Higgs boson. The cross-section for a SM Higgs boson plus 
one jet with Mh = 120 GeV/c^ at lowest order in QCD is approximately 0.1 pb. For tan/3 ~ 30 the cross-section for 
a 120 GeV/c^ pseudoscalar Higgs in the MSSM is about twice as large as for a 120 GeV/c^ SM Higgs at the Tevatron 



8 



10" 



10-2 - 



10-^ - 



10'^ - 



10-^" - 



■|Q-12 I I I I I I I I I I 

50 100 150 200 250 300 350 400 450 500 

Pt [GeV] 

FIG. 3; Transverse momentum distributions for the SM Higgs boson plus one jet production at the LHC with Mh ~ 120 GeV/c^ 
for the different channels. The curves labeled 'Top' ('Bot') include only the top (bottom) quark loops. 

and continues to grow with tan p. The differential cross-section for pseudoscalar plus jet production at the Tevatron 
is shown in Fig. \^ At the Tevatron, the large tan /3 region is completely dominated by bottom quark loops where the 
HQET is of little use. 

For the LHC, the entire tan/3 region is experimentally accessible. In the small tan/3 region, the cross-section is 
well approximated by the HQET limit and the bottom quark contribution can be neglected. However, there are 
regions where both the top and bottom quark loops are important. The results are summarized in Figs. Eland {7\ 
These plots use the full theory matrix elements. For pseudoscalar plus jet production, including only the top quark 
loop underestimates the total cross-section by 9.5% at tan/3 = 4 and the discrepancy becomes larger as tan/3 grows. 
Including only the bottom quark underestimates the total cross-section by 5.6% at tan/3 = 8 and becomes a better 
approximation as tan /3 increases. The total cross-section for the MSSM lightest scalar plus jet production receives an 
important contribution from the interference between the top- and bottom-quark loops over a large range of tan /3. 

VI. CONCLUSIONS 

We calculated the differential distributions and cross-sections for the SM Higgs, H, the MSSM scalar Higgs boson, 
and pseudoscalar boson, A°, plus one jet production at the Tevatron and LHC. We included both the top and 
bottom quark loops and investigated the validity of the Heavy Quark Effective Theory (HQET) limit and the light 
quark mass limit. For large tan/3, the HQET fails and the complete result with all mass dependences is needed. 

The NLO QCD corrections for Higgs plus jet|2^|2^|23| and pseudoscalar plus jet|23 production have been previ- 
ously found in the large mtop limit. Our results make it clear that these can only be applied to the MSSM in certain 
regions. At large tan/3, using the bottom-quark only is a very good approximation in the MSSM. At small tan/3 the 
MSSM pseudoscalar is top-quark loop dominated, whereas the lightest scalar in the MSSM still receives important 
contributions from both the top- and bottom-quarks over a much broader range of tan/3. This can be seen as the 
effective suppression of the cj* coupling and enhancement of the c'^ coupling at small tan/3 where the interference 
between the two terms is still playing an important role. 
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APPENDIX A: COMPLETE PSEUDOSCALAR MATRIX ELEMENTS 



For the qq gA*^ sub-process, the (spin and color averaged) matrix elements squared are particularly simple 
because the presence of a 75 makes the traces much smaller than in the scalar case. They can be written in terms of 
the integrals presented in|2^. 

where the new variables are defined 

si = s-Q^ £i=£-Q^ ui=u-Q'^. (A2) 

It is easy to see that Si — —{i + u) and so on. 

In these expressions we use the notation of jl^]. The Ci loop integral that appears in the calculation is the usual 
triangle integral with two massive legs. For = 0, p^ = = ra|,, pi2 = P1+P2 and = s, the triangle integral is 
defined as 

Ci{s) = Ci{pi,p2) (A3) 
' ^ (A4) 



Z7r2 J [g2 _ rnj] [{q + Pi)"^ - mj] [{q + ^12)^ - 
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FIG. 5: The transverse momentum distributions for the MSSM pseudoscalar Higgs boson for tan 13 = 30, 40, 50 and M^o = 
120 GeV/c? at the Tevatron including the top and bottom quark loops. The top, middle, and bottom lines in the top graph 
are the curves for tan/3 = 50,40,30 respectively. Below is the fraction of the process that comes from each of the different 
channels. The curves at pt = 250 GeV/c from top to bottom are the qg, gg, and qq channels respectively. 



The box integrals with — p'^ = = 0, and 23 = (Pi + P2 + Ps)^ = are defined as 

D{.s,i)^D{p,,p2,P3) (A5) 
' ^ (A6) 



Z7r2 J [g2 _ TO^][(g + - mj][{q + ^12)^ - m'j][{q + pi23)^ - mj] ' 

It is easy to see that the box integrals satisfy the relation D{x, y) = D{y, x). The computer package FF'J^I was used 
to evaluate the scalar integrals. 

For the gg — s- gA'^ sub-process, the (spin and color averaged) matrix element squared can be written in the symmetric 
form. 
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FIG. 6: Cross-section for the production of the MSSM pseudoscalar Higgs boson plus one jet for different values of tan /3 at the 
LHC for M^o = 120 GeV/c^ integrated for pt > 30 GeV/c using the full theory matrix elements. The top and bottom labels 
show what the contribution of the top and bottom quark would be alone. In the region 4 < tan (3 < 8 the total cross-section 
is not represented well by either the top or bottom matrix elements alone. In the experimentally accessible region, the total 
cross-section at the Tevatron is dominated by the bottom loop so only the bottom contribution is shown for tan /? > 30. 
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APPENDIX B: ANALYTIC LIMITS OF MATRIX ELEMENTS 



The partonic cross-section for qq g<^ is 



di 167rs2 36' 
where the spin and color average is explicitly given, 
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FIG. 7: Cross-section for the production of the MSSM lightest scalar Higgs boson plus one jet for different values of tan P 
integrated for pt > 30 GeV/c using the full theory matrix elements. The top and bottom labels show what the contributions of 
the top and bottom quark would be alone. In the experimentally accessible region, the total cross-section at the Tevatron is 
dominated by the bottom loop, so only the bottom contribution is shown for tan/3 > 30. Below is the corresponding mass of 
the lightest scalar for M^o = 120 GeV/c^ . 
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where rrij is the fermion mass in the loop. The integrals are defined by: 



Ilia) — [ dxlogl 1 — ax(l — xj) , /2(a) = / — log [ 1 — aa::(l — a;) 
^ V / Jo X \ 



In the large fermion mass limit, nij — > oo.|21ll2 
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In the small fermion mass limit, to, 0.|22 



11 11 



^qq 



^11 



2m'j 
TOjTr 



qq I 
1 + A J - 



4m^ 



A.-21og(^^ 



1 



A. 



(B7) 



where Kg = log(s/A/^). 

The result for qg q$ can be found from crossing, 
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In the large fermion mass limit, rrij —> oo|2lL l22l|. 
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In the small fermion mass limit, to, — s- 0, 
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(Bll) 



where A„ = \og{\u\l Mfj). 

The results for pseudoscalar production are found assuming the 11:11! A^^ coupling given in Eq. 13 The form factor 
for 



with all moment outgoing and p\ = 0, pi^ — Af^o, [pi + ps)^ = s, is given by 



(B12) 



2t: V 



The differential cross-section for qq gA^ is given by Eq. IBll with 
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In the large fermion mass limit, rrij oo|29j. 
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In the small fermion mass limit, to, — > 0, 
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where = \og{s/M\o). 
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